Photon Aided and Inhibited Tunneling of Photons 
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In the light of the interest in the transport of single photons in arrays of waveguides, fiber couplers, 
photonic crystals, etc., we consider the quantum mechanical process of the tunneling of photons 
through evanescently or otherwise coupled structures. We specifically examine the issue of tunneling 
between two structures when one structure already contains few photons. We demonstrate the 
possibility of both photon aided and inhibited tunneling of photons. The Bosonic nature of photons 
enhances the tunneling probability. We also show how the multiphoton tunneling probability can 
be either enhanced or inhibited due to the presence of photons. We find similar results for the 
higher order tunneling. Finally, we show that the presence of a squeezed field changes the nature of 
tunneling considerably. 

PACS numbers: 42.82.Et, 42.50.-p, 05.60.Gg 



Quantum transport property of Fermions has been ex- 
tensively studied in the last century |l, 2]. One interest- 
ing phenomenon of quantum transport of Fermions is the 
Coulomb blockade effect [3j, which happens when elec- 
trons transport through a quantum dot. The Coulomb 
blockade effect inhibits tunneling by the presence of elec- 
trons in the quantum dot. Because of the Pauli exclusion 
principle for Fermions, each state of electrons my only 
contain one electron; incoming electron must occupy a 
higher energy level. The differences between energy lev- 
els in small quantum dots are very large (with the higher 
level far from the Fermi energy); therefore, the transport 
probability is very small and the tunneling is blocked. 

An analog of the Coulomb blockade [4j, |5( was demon- 
strated by Birnbaum et. al [6] , showing that if an atom 
resonant with a strongly coupled single mode cavity could 
absorb one photon, then the absorption of second pho- 
ton was inhibited. A similar experiment was reported 
in the context of superconducting qubits [7J. Another 
phenomenon is the dipole blockade [8] in Rydberg atoms 
where the excitation of a second atom to a Rydberg State 
is forbidden if one atom is already excited to a Rydberg 
state. The first excitation makes the second excitation 
non-resonant which leads to the blockade effect, and is 
similar to the Coulomb blockade in that it is based on 
the energy gap. The inhibited tunneling of photons that 
we discuss in this paper has a different origin. It not 
only depends on the energy gap but also on the Bosonic 
nature of photons. More generally in the context of trans- 
port of single photons (e.g. through arrays of waveguides 
or through lattices including photonic crystals) one may 
ask how the transport or tunneling characteristics de- 
pend on the presence of photons. Further one may ask 
how tunneling could depend on the quantum statistics of 
photons. 

The simple and clear way to find the answer is to ex- 
amine the photon tunneling between two or more coupled 
modes, with each mode referring to a different structure. 
This can be realized by coupled single-mode waveguide 



devices [9|, [lOj or by fiber couplers. The waveguides cou- 
pled by evanescent fields can be tailored by changing the 
distance between waveguides [9|,[l0( . The evanescent cou- 
pling is responsible for the tunneling of photons. Thus 
the transport of photons in this system is purely deter- 
mined by how photons tunnel from one waveguide to the 
other. 

In this letter we report a remarkable property of the 
tunneling of photons, i.e. the photon aided tunneling 
(PAT) and the photon inhibited tunneling (PIT), which 
can occur due to the presence of photons in the other 
waveguide. We find that even when the energy gap be- 
tween the two waveguides is large, and the single photon 
tunneling is negligible, the PAT is significant and be- 
comes about 1/e when the number of photons in the other 
waveguide is large. The PIT occurs when the energy gap 
is small. When the energy gap is zero, the tunneling rate 
without any photons in the other waveguide is 100%, the 
tunneling is totally inhibited when photons are present in 
the other waveguide. Both the PAT and PIT depend on 
the Bosonic nature of photons, i.e. each state may have 
more than one photon. It is not allowed for electrons due 
to the Fermi statistics. Our results are exact and go far 
beyond the perturbation theory. 

We discuss the tunneling of photons between two cou- 
pled single- mode waveguides labeled A and B (Fig. Q]), 
which can be realized by silica-on-silicon |9|, |lfjj. Ar- 
rays of waveguides have been extensively studied with 
both classical and quantum light |9rll5j. We can also 
use fiber couplers. We start with the simplest possibil- 
ity that the waveguide A contains one photon and the 
waveguide B contains n photons; we calculate the prob- 
ability -P(i.n)->(o,n+i) of the one photon from the waveg- 
uide A tunneling to the waveguide B. If P, h,n)->(o,n+i) > 
-P(i.o)-s-(o.i)> then we conclude that the presence of n pho- 
tons at waveguide B enhances tunneling. Experimen- 
tally, using an APD to detect the probability of having 
zero photons in waveguide A gives us a measurement of 
the tunneling of a single photon to the waveguide B. 




waveguide A 
waveguide B 



FIG. 1: Silica-on-silicon coupled waveguides |g]. 



Note that single photon technology is reasonably well 
developed and we now have an efficient source of single 
heralding photons [16]. Theoretically, the Hamiltonian 
for a system of two coupled waveguides is given by 



H = A (a f a - tfb) + J (a f & + tfc 



(1) 



here a' (tf) is the operator that denotes creating a pho- 
ton in waveguide A (waveguide B), and a (b) is the 
corresponding annihilation operator. Generally, the two 
waveguides A and B have different refractive indices 
Eq ± A. In the discussion of the time evolution of states, 
£o only gives an overall phase e _le °', which does not af- 
fect the transport probability. Thus, without loss of gen- 
erality, we set £o = 0. The coupling constant J is the 
tunneling energy of the photon to go from one waveguide 
to other. 

If the system is initially in the state <I> (0) = \n, m) , i.e. 
there are n photons in the waveguide A, and m photons 
in the waveguide B, then at time t the system would be 
in state 



* (t) = V JJ\_K _J7_ | 0) Q) 

\/n\m\ 



(2) 



The time dependent operators a'{— t) and b*(—t) can 
be obtained through the Heisenberg equation of mo- 



tions. We have 



ot(-t) 

bH-t) 



V(t) 



(0) 



where 



V{t)=e 



-iMt 



with M = 



&t( ) 

The V (t) can be 



A J 

J -A 
directly calculated from e 1CT ' A — ±i wo^i T n,i/ ■ **■) — p 
Here a is the pauli matrix, A is the magnitude of vector 
A, I2 is the two dimensional identity matrix. We can 
then find 



I 2 cosA + y{cj-A) s ^. 



V(t) 



-k/P 



-iVP 



(3) 



On defining 7 = — , Q = \J\ +j 2 and Pq 

amplitude of the off-diagonal term of V is given by 
with 

P = P sin 2 {QJt) . 



■§2, the 



The amplitude of diagonal term is given by \/l — P, and 
the corresponding phase 6 can be calculated from prob- 
ability yT^Pe- ie = cos (QJt) - i^ sin (QJt). 
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FIG. 2: The tunneling probabilities -P(i,o)->(o,i) (red dashed) 
and -P(i,i)_).(o,2) (blue) as a function of (a). P; (b)-(d) scaled 
time Q^L for different values of 7 = 4 • 



At time t, the probability that the system is in the state 
\n',m') is given by P( n , m )^y{n> ,m<) = \(n' , m'|* (i))| 2 . If 
the system is initially in the state V I / (0) = |l,n>), i.e. 
one photon in the waveguide A and n photons in the 
waveguide B, then the tunneling probability of the one 
photon from A to B is given by 



P(l,n)^(0,n + 1) = (ft + 1) (1 - P) U P 



(5) 



The standard case of tunneling is a special case of 
Eq.([5]) when there is no photon in waveguide B, i.e. 
n = 0. We may have -P(i,o)-s.(o,i) — P> (For brevity, 
we use P instead of -P(i,o)-s-(o,i) m the later discussion). 
Notice that this one-particle tunneling probability P is 
the same as in the Fermion case. In fact, we can see that, 
when A is large compared to J, then the energy gap be- 
tween two modes is large (which means a high barrier at 
the junction), leading to large 7 and a small tunneling 
probability. The small tunneling between the large gap is 
in fact the theoretical precondition of the Coulomb block- 
ade of Fermions, it does not allow Fermions to tunnel to 
a state with much higher/lower than the Fermi energy. 

In order to illustrate how the presence of photons 
within the other waveguide affects tunneling process, we 



fp first compare the tunneling probability P(in^.(o,2) with 
the presence of one photon in the other waveguide with 
P. From Fig. [2] (a), we can see that when P < 1/2, 

(4) -P(i,i)^(o,2) > P, i- c - the presence of one photon in the 
waveguide B enhances the probability of tunneling, i.e. 
PAT occurs. It is easy to show from Eq.([5]) that the max- 
imum difference (-P(i,i)->(o,2) — P)mnx = 1/9 is reached 
when P = 1/3. This means that the PAT occurs for any 



range of values of the A and J. When the gap 2 A is such 
that 7 = y > 1, we have P < 1/2, P < 1/2, then the 
PAT always occurs. This is shown in Fig. [21(b). Here 
7 = 2, Pn 1)^(0.2) > P f° r an Y time t. When the energy 
gap is small 7 < 1 thus Pq > 1/2, we can still observe 
the PAT. The oscillation structure of P, Eq.Q guaran- 
tees that there exist time regions such that P < 1/2 so 
that P(i,i)=-(o,2) > P. (Fig. [2(c) and Fig. [2(d)). 

From Fig. [2(a), we also observe the appearance of the 
photon inhibited tunneling (PIT) when P > 1/2. Espe- 
cially when P = 1 (100% tunneling probability without 
present any photon in waveguide P), we have exactly 
f(i.i)->(o.2) = 0, the photon tunneling is totally inhib- 
ited. The PIT can only be observed when the gap is 
small 7 < 1 so that Pq > 1/2, which allows P > 1/2 in 
some time region (Fig. [2(d))). 

The existence of PAT and PIT shows a competition 
mechanism introduced by the Bosonic nature of photons, 
as is seen from Eq.©. Two factors are multiplied to 
P: the first one (n + 1) > 1 is an aided term, which 
makes P(i,n)->.(o.ri+i) possible to be bigger than P. It 
is from the Bosonic nature that a) \n) = \Jn + 1 \n + 1). 
The Bosonic nature shows that if a state contains more 
Bosons, then it is easier to add an extra Boson in it. The 
second term (1 — P) < 1 is an inhibited term, which 
reflects the tendency that the n photons from waveguide 
B stay in the waveguide B. Notice that (1 — P) gives 
the probability that the one photon still wants to stay 
at its own site. The power n reflects the probability of 
the n-photon state to remain n-photon state. When n 
is larger, this term becomes smaller. The aided term 
(n + 1) is fixed when the photon number is fixed; while 
the inhibited term (1 — P) n can decrease from 1 to 
when P increases from to 1. As a result: if P is small, 
the aided terms dominates and we observe the PAT; if P 
is large, then we observe the PIT. 

In the light of this discussion, we examine P(i,i)->.(o,2)- 
When P — 0, every photon must stay in its own state, the 
tunneling probability should be zero. When P < 1/2, the 
tunneling probability is small, every photon still wants to 
stay in its own state, the inhibited term (1 — P) n is large. 
At this time the aided term (n+ 1) is important, leading 
to -P(i. 1)^(0.2) > P- When P > 1/2, the inhibited term 
suppresses the positive effect of the aided term (n + 1). 

We now discuss the general case of the tunneling prob- 
ability -P(i,n)->-(o,n+i) of one photon tunneling in the pres- 
ence of n photons. From Fig[3](a), we can see that 
when number of photons is increased, the region of PAT 
becomes smaller and occurs at the smaller values of 
P, however the maximum value does not decrease very 
much. In fact, from Eq.©, we get (P(i,„)^.(o,n+i))ma X = 
(1 — ^zj)™ when P = l/(n + 1). When the number of 
photons in waveguide B is very large n — > 00, we can get 
(• p (i,n)^(o,n+i))max -► 1/e ~ 0.37. Therefore we get very 
significant PAT especially when A is large which leads to 
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FIG. 3: (a). The tunneling probability P(i, n )_>.(o,n+i) as a 
function of P for n = (red dashed); n = 5 (green); n = 10 
(black); n = 50 (blue), (b)-(d): The time dependence of the 
tunneling probability -P(i,5o)-+(o,5i) (blue) as compared to that 
of P(i, o)^(o,i) (red dashed) for different values of 7. 



very small values of P. We compare the tunneling prob- 
ability P(i.„)-s.(o,n+i) w hh P for different 7. In Fig[3](b), 
when 7 = sfn = \/50 so that Po = 1/51, we can see 
that the tunneling without any presence of photons in 
the waveguide B is very small, while P(i.5o)=-(o,5i) h as 
a plateau at 0.37 for a large time range; In Fig[3](c), 
the aided tunneling is still significant for 7 = 5. In 
FiglSl(d), when the gap is small, the aiding tunneling 
acts as the pulse in the vicinity of time period T . This 
remarkable tunneling for large n has important appli- 
cation. This means that even for a large gap A, we 
can always find the finite photon tunneling probability 
(•P(i,n)-Ko,«+i))ma X near ± by choosing n > -y 2 = ^j-, 
whereas as for large A, P(i,o)->(o,i) is negligible. 

The above PAT can be generalized to multi-photon 



tunneling. In Figf4l we compare P( 
P" 2 . Obviously, P, 



(712, n)— »(0,n+n2) 



with 



(n 2 ,0)^(0,n 2 ) 



(n 2 ,0)->(0,n 2 ) 



decreases 



when ri2 increases since P < 1. This is intuitive in that it 
is harder for more photons to tunnel to another waveg- 
uide. We can see that when 712 = 10, the probability 



(n 2 ,0)-K0,7i 2 ) 



is close to zero for most values of P. How- 



ever, with photons in the waveguide B, the tunneling 
is significant. We can always find a finite maximum of 
P(„ 2i „)_i.(o,n+n 2 )- If n 3> n 2 , this tunneling probability is 
even greater than the one-photon tunneling P(i j o)=.(o.i) 
as shown in Figf4j 

Using Eq. ([2]), we calculate the result for 



P(n 2 ,n)->(0,n+n 2 ) to be: 



P 



(n 2 ,n)— >(0,n+n 2 ) 



n 2 



1)2 



(1-P)"P" 



(6) 
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P{n-2,n)^,(o,n+n 2 ) (thick lines) in presence of 10 photons 
in the waveguide B, as a function of P for 112 = 2, 3, 5, 10. 
For each line, ni is marked at it's typical peak. The multi- 
photon tunneling probabilities without presence of photons 



in the waveguide B, P ( 
lines. 



where 



(n 2 ,0)->(0,n 2 ) 



P" 2 are given by thin 



t — ^V; — r is the binomial coefficient. 

m I (n— m)!ml 

From this equation, we can clearly see that the inhib- 
ited term is still given by (1 — P)™, while the aided 

n + n 2 

n 2 
tunneling should be more important when n 2 increases. 



term is changed to 



This means the aided 



However, noticing that P( 



n 2 ,0)-K0,rc 2 ) 



decreases, the ab- 



solute value -P(n 2 ,n)-j-(o,n+n 2 ) mav still decrease. This can 
be seen from Fig|4] which shows the decreasing peak of 



(n 2 ,n) — >(0,n+n 2 ) 



with the increase of n 2 . The maximum 



peak value is given by 



P 



(n 2 ,n)— y(0,n+n 2 )j 



n + n 2 
n 2 



(n + n 2 ) 



n+n 2 



(7) 



which is reached when P = 



" 2 



The limit when 

n+n 2 

is the maximum of -J-i^^™ 2 



n — >• oo, is tiie maximum ol -^-l -)" ! . It decreases 

' n 2 \ v c ' 

when n 2 increases. However, even for n 2 = 10, we 
still have the finite tunneling about 12.5%, which is 
much greater then the corresponding tunneling proba- 
bility -P(io,o)->(o,io) — (^o) without the presence of any 
photons in the waveguide B. 

In an experiment, it is much easier to prepare the field 
in a coherent state than a state with fixed photon num- 
ber. The coherent state |/3) = e^l 2 / 2 XT=o y=T I") has 
the average photon number n — (f3\tfb\/3) — \(3\ 2 . We 
may then discuss the possibility of the PAT and PIT 
with the field in coherent state in the waveguide B with 
a fixed average photon number h. The probability of 
P(„ 2;( g)_j.(o ; /3,n 2 ): n 2 photons tunneling from waveguide A 
to the waveguide B is given in terms of the confluent 
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FIG. 5: Comparison of photon tunneling probabilities with 
field within waveguide B in a coherent state and in a 
Fock state, (a). The single photon tunneling probability 
P(i ;/ 3)^(o ; /3,i) (solid lines) and P(i ; „)-K0;n+i) (dashed lines) 
as a function of P for n, n = 1,10; (b). The multipho- 
ton tunneling probability P(„ 2 ;/3)^,(o;P,n 2 ) (solid lines) and 
P(n 2 ,n)->(o,n+n 2 ) (dashed lines) in presence of (average) 10 
photons in waveguide B, as a function of P, for n-2 — 2, 10. 
The black dashed dotted lines represent tunneling probability 
for a field in squeezed vacuum in waveguide B. 



hypergeometric function \F\ [a, b;z\. 



(ra 2 ;/3)-K0;/3,n 2 



1F1 [l + n 2 ,l,n(l-P)]P" 



(8) 

In the special case of the one photon tunneling, the result 
is rather simple 



■P(i ; /3)->(0;/3,i) 



-nP 



[1 + n (1 - P)] P. 



(9) 



The formula for P(itf)-+(o;p,i) and P(„ 2;/ s)-K0;/3,n 2 ) are 
complicated, it's not easy to separate aided and inhib- 
ited terms. It's no doubt that e~ n and e~ nP are inhibited 
terms as they are always smaller than 1 ; however the sec- 
ond term contains both aided and inhibited parts. Even 
for one photon tunneling, the factor [1 + n (1 — P)] can 
be either greater or smaller than 1. However, from FigJSJ 
we observe that, the results for the coherent state case are 
very similar to the case of a Fock state in waveguide B, es- 
pecially when n — n is large. We examine the maximum 



ofP (1 . 



/3)-HO;/3,l)> 



which is 



Vft 2 +2n+5-2 ^ (y / n 2 +2n+5-ft-3) 



and OCCUrS for P = " +3ft-nVn 2 + 2ft+5 



For large n, 

linifi^oo (P(i ; /3)- + (o:/?,i)) max = ~> which is same as for the 
Fock state. 

The tunneling probabilities are also sensitive to the 
photon statistics of photons in the waveguide B. To il- 
lustrate this, we consider the field in a squeezed state 
l 



10 



Vc 



Eoo 
»=0 



,imp 



(tanhr) 



n!2™ 



|2n) with n = 



sinh r = 10. In this case, calculations show that 



P 



(n 2 ;£)-K0;£,n 2 ) 

pn2 .*[! 



(10) 



n 2 2 



VT 



n 



n 2 , „ N 9 n 
2 1 ; (1_P) 2 



where 2 F\ [a, b, c; z] is the hypergeometric function. The 
dash-dotted black line in the Fig JS] shows the behavior of 



(110)) for ri2 = 1 and 10. The behavior is clearly different 
from the case of a coherent state: a long plateau occurs. 
We can conclude that the tunneling probability in the 
presence of a field in squeezed state -P(n 2 ;?)-K0;£,n 2 ) i s 
mostly inhibited compared to the case when no field is 
present in waveguide B. 

In conclusion we have shown how the tunneling of a 
single photon as well as multiphoton tunneling can be 
enhanced or inhibited by the presence of photons. We 
presented the physical reasons behind such an enhance- 
ment or inhibition. A crucial role is played by the Bosonic 
nature of photons. The waveguide structures or fiber cou- 
plers are known to be almost decoherence free, however 
if need be then the decoherence effects can be taken into 
account using the formulation of Ref . [17( . Although we 
explicitly considered the simplest case of a coupler the 
results can be extended to arrays of couplers. Tunneling 
is a universal effect in Physics; and therefore, results of 
this letter would be applicable to all situations in which 
Bosons are involved. Further the results of this paper 
should have a bearing on the quantum walk of a single 
photon in presence of other photons [181 ] . 

Xuele Liu would like to acknowledge Amanda Taylor 
for a careful reading of the manuscript. 
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